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1 Introduction 



The phenomenology of b hadrons is very rich, because owing to the large mass of the bottom quark there 
are many decay channels. Available semileptonic decays into any hadronic part containing a charm quark 
are among the decay modes of B mesons which play a critical role in the analysis of unitarity triangle 
and are an essential ingredient to any precise information about the CKM matrix element \V c b\. 
However, there are many puzzling features associated with the semileptonic b — > c data which have existed 
during the last ten years as, for example, the so called "1/2 versus 3/2 puzzle" [1,2]. It corresponds to 
the persistent conflict between theoretical predictions and experimental measurements of semileptonic 
branching ratios of B — > D**iv: 

T(B d ->■ Dl* /2b Jv) < T{B d -> D* 3 * /2 iv) [Theory] 

T(B d ->■ D\* /2 iv) > T(B d -> D'** /2 Iv) [Experimental data] 

where D** are the first orbital excitations of D mesons having a positive parity. In the quark model 
picture, if we consider P wave states we see that four orbital excitations of heavy mesons should exist. 
In the case where the heavy quark is the charm quark c, we obtain the D** represented in table 1: So, 
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Table 1: Classification of D** states. The total angular momentum J = j '■ + s c is the sum of j (the angular momentum of 
the light component of a heavy meson) and s c = 1/2 (spin of the heavy quark). The first column contains the classification in 
the infinite mass limit where due to the heavy quark symmetries, j is a conserved quantity and thus a good quantum number to 

classify the states. 



the comparison between experimental and theoretical results presents some hardships and this challenges 
our understanding of QCD and is of high relevance for an accurate estimate of V c b- 

However, in many theoretical approaches (HQET, heavy quark expansion, quark model, Lattice QCD 
with quenched approximation, etc.), branching ratios corresponding to the B — > D** iv decay were 
calculated using the infinite mass limit. That is the reason why, in order to address the aforementioned 
questions, we propose to determine the form factors and then the branching ratios using "real" charmed 
quarks having a finite mass. 

In the following, I will explain the extraction of form factors corresponding to the decay of B — > D** in 
the framework of Lattice QCD and I will discuss some preliminary results we have obtained. This work 
is done in collaboration with B. Blossier and O. Pene (LPT Orsay France), V. Morenas (LPC Clermont 
Ferrand France) and K. Petrov (LAL Orsay France). 

2 Form factors 



We choose to consider in our work the decay of B into the scalar 3 Pq and the tensor 3 P2 states of D** . 
The semileptonic decay of a pseudoscalar meson into a scalar meson is mediated by the axial part of the 
weak (V — A)^ current. The operator denotes the vector current qj^b and the operator A^ represents 
the axial current 97^756. 

The matrix element for the 3 Pq state can be parametrized in terms of two form factors (u+, 



( 3 Po(p D *,) I A,* j B(p B )) = u+ {p B +p D „)n + u- (p B -p D „) 



and for 3 P2 state: (Here, A is the polarization tensor corresponding to J = 2 states, A = {0,±1,±2}) 



( 3 P 2 fe„, A) I F M I B(p B )) = 
( 3 P 2 (p D „ 1 \)\A,\B(p B ))=\I\ 
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3 Going to the Lattice 



The determination of many observables such as form factors, decay constants as well as numerous matrix 
elements, which play an important role in Flavor Physics, requires the use of non perturbative methods 
because when we look at the low properties of QCD, we can no longer use perturbative techniques 1 . 
Lattice QCD 2 is considered as the only way to systematically and rigorously solve non perturbatively 
the quantum theory of strong interactions starting from first principles. 

Briefly, it is a means of regularizing Field Theory where the continuum and infinite space-time is replaced 
with a discretized grid of points in a finite volume of extent L in space and T in time, separated by a 
distance a (i.e. the lattice spacing). Quark fields are present on the sites of the lattice and gauge fields 
are the links between those sites. 

3.1 LQCD action 

The gauge action used in our simulation is tree-level Symanzik improved [6] with j3 = 3.9 corresponding 
to a lattice spacing a — 0.0855 fm and where b\ = —1/12, bo = 1 — 8bi. 

The fermionic action is the Wilson Twisted-mass Lattice QCD (tmLQCD) action with two flavors of 
mass-degenerate quarks, tuned at maximal twist in the way described in full details in Ref. [5]: 



3.2 Simulation setup 

We have computed green functions using two ensembles of gauge configurations produced by the European 
Twisted Mass Collaboration (ETMC). Simulation parameters are presented in Table 2. 






L 3 x T 




Mc 




nb. of gauge configurations 


3.90 
4.05 


24 3 x 48 
32 3 x 64 


0.0085 
0.0060 


0.215 
0.1849 


0.3498, 0.4839, 0.6694 
0.3008, 0.4162, 0.57757 


240 
200 



Table 2: Parameters we have considered in our simulation. ft,/i C i|Uft are respectively the bare twisted light, charm and heavy 

quark masses 



4 Determination of form factors 

The first step in the determination of form factors, presented in section 2, is the computation of matrix 
elements contributing to these form factors, and in order to access matrix elements on the lattice one 
computes the following three-point correlation functions 
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Figure 1: Sketch of the valence quark flow in the form factor of B — > D**£v 

^{t,t u t f ,p u p f ) = J2 (&D*.(tf,$f) M*>Z) ^B(U,Xi)) ■ e *< 2 -*/>-fc • e-*( 2 -*)-ft 

positions 

where vL**, &b are respectively the creation and annihilation operators of D** and B mesons, is the 
vector or axial current. 

*At short distances or at high energies the quarks interact weakly, so that it is possible to study the theory of strong 
interactions (QCD) with perturbative techniques since the coupling constant Qq CD <C 1 
2 For further reading, I refer the reader to [3,4] 



From the asymptotic behavior of the three-point correlation function, it is clear that the removal of the 
exponential factors can be achieved by considering the ratio 

m = m y(3) ^'/ffft> VS5-V35 (1) 

where S'm — |(0|^m|-^0| 2 is obtained from the fit with asymptotic behavior of the two-point correlation 
functions. 

tfC0( t ) > e -EMt 

t— >oo 2E]\j 

When the operators in the ratio (1) are sufficiently separated in time, one observes the stable signal 
(plateau), which is the desired hadronic matrix element: 

4-1 k form factor 

In order to obtain simpler expressions representing the interpolating fields ! of the P wave excited states, 
we choose to work in the rest frame of D** mesons where p D ,, = (m c „ , 0) and we also choose a symmetric 
momentum for the B meson p B = (E B ,p,p,p). Using different values of the recoil parameter w = € 

{1, 1.025, 1.05, 1.1, 1.15, 1.2, 1.3}, we perform the extraction of k which is the form factor contributing 
the most to the decay of B into D** . 

Combining the three point correlation functions contributing to k 4 , we study the difference with respect 
to the infinite mass limit ones: 

(3) ~ ] ) 

^infinite mass limit = P k x ' == 

where k x is the form factor at the infinite mass limit. Using what was found in [8]: 

i = V3 7^7 (1+w) t 3/2 (w) m D , = r D , m B 

2 a 2 

r 3/ »=r 3/2 (l) 3/2 r 3/2 (l)~ 0.539 and ^-1.50 

we find that for w — 1.3 (using for the moment (3 = 3.90): 



7 finite 



c*(3) 

^infinit 



(3) [ afi = 0.34 ~> 1.38 ±0.49 

a^ = 0.45 ~* 1.97 ±0.65 
an = 0.67 ~+ 2.88 ± 3.82 



The correlators corresponding to the higher b mass are characterized by large fluctuations in the effective 
mass plateaus and thus by large statistical uncertainties. Taking these results with a pinch of salt, and 
excluding data of the highest heavy quark mass from the next analysis, we estimate roughly the branching 
fractions MB -> D , ) ■ MB , -> D*ir). 

\ (2+)' y (2+) ' 

SUB - D + ) .*(D + - DV) = = 2 ' 4 ^ (3 - 6 ± 2 - 6) X 10 " 

V < 2+ ' ; V < 2+ » ' |m BGoV =2.9 ~* (7.4 ±4.9) x 10~ 3 

Although we are still far from the physical B mass, we have an indication that there are many observables 
in B physics that can be determined using LQCD. When results are more refined we aim to do a more 
detailed comparison with the experimental data. Results will be published as soon as ready. 



3 interpolating fields are the mesonic creation and annihilation operators which defines the quantum number of a state. 
To find the interpolating fields of 2+ states, we followed a strategy based on group theory where non local operators are 
implemented [13]. 

4 For example, if we consider the A = polarization of the 3 P2 state, k could be written as: 



-f( 3 P 2 (0)\A 1 \B(p B )) 



4-2 Scalar transition amplitude 

The determination of matrix elements corresponding to the decay of B into 3 Po(p D »» ) state, at different 
recoils w, leads to the weak form factors u + and u_. For the moment, we focus our attention on the 
extraction of ( 3 Po(p D ,, ) | Aq | B(p B )) at zero recoil. It seems that this hadronic matrix element is not 
equal to zero contrary to what was already found in the infinite mass limit. Note that it is an important 
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Figure 2: ( 3 P (m D „ ) | A \ B(m B )) at = 3.9 



Figure 3: ( 3 P (m D « ) | A \ B(m B )) at @ = 4.05 



result that the matrix element ( 3 Po(m D „ ) | Aq | B(m B )) is non vanishing. It shows us how working with 
a realistic charm quark could lead to much more accurate results than working with the infinite mass 
limit approximation. 



4-3 ^o(l) a t zero recoil 

As a byproduct of our analysis, we can also study the B — > D* Iv decay at zero recoil. The corresponding 
form factor reads [11,12]: 

where is the axial renormalisation constant taken from [9]. The values obtained for J^o(l) are presented 
in tables 3 and 4. They do not show any contradiction with what was already found using Lattice QCD 
with a different fermionic action and different approximations, but the extrapolated value (Fig. 4) seems 
to be higher than what was already found using QCD sum rules [10]. 



a/it 


^o(l) 


0.35 


0.882(44) 


0.48 


0.874(48) 


0.67 


0.888(58) 



afj,b 


^o(l) 


0.301 


0.814(48) 


0.416 


0.821(50) 


0.5757 


0.827(53) 



Table 3: ^b(l) at /3 = 3.9 for different values of a[t b 



Table 4: ^o(l) at /3 = 4.05 for different values of a^ b 



5 Discussion 

During the last ten years, techniques and algorithms in LQCD have been developed and computing power 
has increased: this is really encouraging for the Lattice community where some important parameters in 
B physics are and will be calculated reaching a precision on par with the experimental one. 
However, working with excited states on the lattice is not a trivial issue. On one side, analytical expres- 
sions of B — > D** transition amplitudes, form factors and decay rates can be calculated. On the other 
side, isolation of such excited states is very delicate and there is an increase in noise when going to high 
momentum and discretization errors become larger when the b quark mass increases. We hope that with 
higher statistics, results will be more clear in order to extract the sought form factors and to discuss their 
implications. 



Figure 4: Extrapolation of J?o(l) at finite physical b quark masses 
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